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The effects of thermal fluctuations of the mass (horizon area) and electric charge, on the entropy 
of non-rotating charged macroscopic black holes, are analyzed using a grand canonical ensemble. 
Restricting to Gaussian fluctuations around equilibrium, and assuming a power law type of relation 
between the black hole mass, charge and horizon area, characterized by two real positive indices, 
the grand canonical entropy is shown to acquire a logarithmic correction with a positive coefficient 
proportional to the sum of the indices. However, the root mean squared fluctuations of mass and 
charge relative to the mean values of these quantities turn out to be independent of the details of the 
assumed mass-area relation. We also comment on possible cancellation between log (area) corrections 
arising due to fixed area quantum spacetime fluctuations and that due to thermal fluctuations of the 
area and other quantities. 
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The notion of entropy of a spacetime such as a black hole presumes the existence of microstates whose degeneracy g 
XJ^ 1 is measured by that entropy. It is conceptually difficult to think of microstates of a classical stationary spacetime. Yet 
| even vacuum solutions of Einstein's equation like the Schwarzschild and the Kerr black holes are supposed to possess 
■ an entropy equal to a quarter of their horizon area (in Planckian units) [1], [2], The only way in which microstates 
might arise in these cases is within a quantum mechanical description of such spacetimes themselves (and not merely 
of the matter fields propagating in them) which allow for dynamical fluctuations. In other words, the entropy of black 
^ ' holes must be taken to originate from quantum degeneracy of spacetime geometry, and consequently, can only be 
r*"** , properly analyzed within a quantum theory of gravitation [1] . 

Of all attempts to construct a 'quantum general relativity', one that perhaps comes closest to dealing directly with 
fluctuations of spacetime geometry (within a canonical framework) is Loop Quantum Gravity (LQG) [3]. Yet even 
in this approach, the only situation that can be reliably analyzed is one in which the black hole is isolated, i.e., has 
a horizon with a fixed (macroscopic) area A. » 1 (in units of Planck area). From a statistical thermodynamics 
standpoint, this corresponds to employing a microcanonical ensemble to compute the entropy [4]. Microstates leading 
q | to the microcanonical entropy Smc = log g are identified with boundary states of a three dimensional Chern-Simons 
theory 'living' on the horizon [4]. The local Lorentz invariance is gauge-fixed to a local SU (2) on an appropriate spatial 
slice. Thus one may take this SU(2) to be the gauge group of the Chern Simons theory. However, boundary conditions 
(3J[), associated with the 'isolatedness' of the horizon have been argued [4] to further gauge-fix the SU(2) to a residual 
U(l). In either case, the degeneracy of the boundary states of the Chern Simons theory yields (for macroscopically 
large horizon areas) a term logarithmic in the area together with an infinite series in inverse powers of the area [5] , 
[6], over and above the Bekenstein-Hawking entropy Sbh — A/ '4 (in Planckian units), 



Smc = Sbh -\\ogS BH + const. + 0{S B ]j) for 5(7(2) 

= S BH -\ log S BH + const. + 0(S BH ) for U(l) (1) 

The negative sign of the dominant logarithmic correction in both cases is due to the fact that invariance of the horizon 
states under the gauge group effectively reduces the degeneracy as compared to the leading order term (area law) 
where the invariance is not taken ito account. Notice that each term of the infinite series in eq. (1) is finite and 
unambiguously calculable. 

The problem with an isolated black hole is that it neither radiates nor absorbs anything. This is unsatisfactory from 
a physical standpoint since classically a black hole absorbs everything, and with quantum matter around, it radiates 
all kinds of quanta in a Planckian distribution [2] at the Hawking temperature given by the surface gravity on the 
horizon. Thus, the appropriate equilibrium ensemble to probe black hole thermodyamics ought to be the canonical 
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(Gibbs) ensemble (generalized to the grand canonical ensemble when fluctuations of charge or angular momentum 
are included). Unfortunately, however, for asymptotically flat black hole spacetimes, the equilibrium is unstable, 
because of a superexponential growth of the density of states leading to a diverging canonical partition function. So, 
depending on whether the ambient temperature is larger or smaller than the Hawking temperature, the black hole 
either accretes or radiates without limit. 

The situation is more favourable for asymptotically anti-de Sitter black holes which, for a certain range of pa- 
rameters, can indeed be in stable thermal equilibrium [7]. For such spacetimes the area (energy) fluctuations within 
a canonical ensemble lead, within a saddle-point approximation, to an additional correction term in the canonical 
entropy, over and above the infinite series of corrections to the area law found for the microcanonical entropy [8] , 

SthSc = ^ogS BH - ilog(n-l) + const, (2) 

where n is the index appearing in the assumed power law relation between area of the horizon and energy. This 
additional correction has a positive sign, as expected from the fact that thermal fluctuations tend to disorder the system 
and therefore always increase the canonical entropy. The area-dependent part is universal in that it is independent 
of the index characterizing the asuumed mass-area power-law relation. We may remark that eq. (2) is valid in the 
limit of large horizon areas and for fixed (i.e. thermally non-fluctuating) charge or angular momentum of the black 
hole. For asymptotically flat black holes, n = 1/2, implying that the entropy as well as the free energy acquires an 
imaginary part. This indicates a breakdown of the saddle-point approximation used to compute the canonical entropy 
from the canonical partition function. It is not unlikely that this breakdown is indicative of the thermal instability 
ensuing in the canonical ensemble for such spacetimes, as mentioned earlier. On the other hand, for all anti-de Sitter 
black holes, n = 3/2, and therefore the calculation above is reliable. 

One immediate fallout of this additional correction is that it begins to compete with the microcanonical corrections 
to the area law. The net canonical entropy now takes the form 

Sc = Smc + S t h Sc 

= Sbh - logSW - -log(ra-l) + const. + ... . (3) 

for the SU{2) case, while for the U{1) case, the two leading logarithmic corrections to the area law simply cancel each 
other out [9] . At this point the cancellation is not much more than a point of curiousity, since there is still an infinite 
series of other (power- law) corrections hidden in the ... in eq. (3). However, if such cancellations are proven to work 
for all orders in powers of inverse area, one may ponder about symmetries, non-renormalization theorems etc. 

Recall that so far the discussion of the canonical entropy has assumed fixed charges/angular momenta of the black 
holes. In this paper we consider the effect of including thermal fluctuations of the electric charge for non-rotating 
spherically symmetric black hole spacetimes. In other words, we consider the generalization of the previous canonical 
ensemble formulation to a grand canonical ensemble framework where the electrostatic potential of the black hole plays 
the role of a chemical potential. Thus, the heat bath which the black hole interacts with is now electrically charged. 
The use of a saddle-point approximation in this case involves important technical changes which are incoporated in 
our treatment. The object here is to determine whether the universality found in the canonical treatment survives 
the generalization. Once again, we make no more assumption than is absolutely essential. In particular, we adhere to 
the area spectrum strictly as derived in the LQG formulation [3], specialized to the case of large macroscopic areas 
where the spectrum becomes equally spaced. 1 

In this context, we may note that two recent papers [13], [14] have explored the effect of thermal energy and 
charge fluctuations for a grand canonical ensemble of AdS Reissner-Nordstrom black holes. While our result for the 
leading thermal fluctuation correction to the grand canonical entropy for this particular black hole matches with 
the results quoted in [13], the implied assumption in those papers regarding the independence of energy and charge 
fluctuations, as also about their relative size, are not tenable even for large black holes, as we shall show in what 
follows. Consequently, the conjectured complete cancellation (including angular momentum fluctuations) between the 
microcanonical and thermal logarithmic corrections to the area law, as proposed in [14], appears to be without basis. 

The paper is organized as follows: in section 2 we present a grand canonical ensemble formulation of small fluctu- 
ations for general equilibrium statistical mechanical systems. This is then used in section 3 to compute the leading 



1 Only in one other approach, involving a proposal for coherent states [10] within LQG, is there a proper derivation of an 
equally-spaced area spectrum [11] which does not involve ad hoc assumptions regarding the black hole mass spectrum. Such 
ill-founded assumptions often lead to claims regarding corrections to the area law [12] which we find spurious. 
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logarithmic corrections to the grand canonical entropy, via a Poisson resummation procedure followed by a saddle 
point approximation to perform the integrals in the grand canonical partition function. In section 4. we study the 
relative sizes of fluctuations in more detail to show that they indeed ensure that the approximation of small Gaussian 
fluctuations is valid for large area black holes, when the charge is also allowed to grow with area in a certain way. 
However that is shown not to be true for black holes of vanishingly small charge considered in [13]. We discuss our 
results and conclude in section 5. 

II. MASS AND CHARGE FLUCTUATIONS IN THE GRAND CANONICAL ENSEMBLE 

For a thermodynamic system with Hamiltonian H and charge operator Q in equilibrium with a heat and charge 
reservoir at temperature T and electrostatic potential $, the grand partition function is 

Z G = Tr exp-/3(# - $Q) (4) 

When the spectra of H and Q are both continuous with eigenvalues E and Q respectively, the partition sum in eq. 
(4) can be conveniently rewritten as an integral 

Z G = J dE dQ p(E, Q) exp-0(E - $ Q) , (5) 

where the density of states p(E, Q) may be taken to be a smooth function and (3 = T -1 . 
Defining the microcanonical entropy Smc(E, Q) = log p(E, Q) 2 we get 

Z G = J dE dQ cxp {S MC -(3{E-<S>Q)} , (6) 

Assuming that the integral is dominated by the saddle point of the exponent at E = M, Q = Q 0} we expand the 
exponential in the integrand upto terms quadratic in the fluctuations SE = E — M , SQ = Q — Q n 

Smc(E, Q) ~p(E-$Q) = S MC (M, Qo) - (3(M - $ Q ) 

+ 2 ( s mc,ee) m ^ Qo 8E 2 + {Smc,qq) m ^ Qo SQ 2 + 2 (Smc.eq) m,q SES Q > ( 7 ) 

where, Smc,ee = d 2 SMc/dE 2 , etc. In writing eq. (7) above we have made use of the saddle point conditions 

(Smc,E) MQo = (3 = T _1 

(Smc,Q) m , Qo = -/?$ = -$/T (8) 

which are the usual relations for the temperature and potential in the microcanonical ensemble. 

The presence of cross terms involving SE and SQ in the exponent in eq.(6) (using (7) implies that the energy 
and charge fluctuations cannot, in general, be treated independently. Evaluating the Gaussian integrals over the 
fluctuations yields, 

2ir 

Zg = -7== exp {S M c(M,Qo)-pM + (3<S>Q } , (9) 
Vdet it 

where, the Hessian matrix 

n = ( S Q MC ' EE s c mc ' e q) (io) 

\&MC,EQ b MC,QQj MQa 

is assumed to be negative definite to ensure stability under small fluctuations. The necessary and sufficient conditions 
for this are 



2 Observe that Smc defined as above is precisely what we called Smc in [8] where we defined Smc = log g(E) with g(E) being 
the degeneracy of states with energy E. This distinction between the two definitions of microcanonical entropy is relevant for 
logarithmic corrections to the area law [8] . 
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TrQ, = Smc,ee\m,q + Smc,qq\m,q < ® 

det ft = [Smc.ee Smc.qq ~ s mc,eq] m ,q > ( n ) 

which necessarily imply 

Smc,ee\m,q < 0, and Smc,qq\m,q < 0. (12) 

Note that while these conditions together imply the first of the necessary and sufficient conditions (11) for stability 
of ft, they are not sufficient to guarantee the second one. 

Using the microcanonical relations for temperature and potential, we may expess the necessary conditions for 
stability in terms of the heat capacity Cq = (dE/dT)Q and the capacitance C = (<9Q/<9$)_e in the following way 

C Q > and C-f < T (^) . (13) 



E 



The more stringent necessary and sufficient conditions can also be similarly expressed in terms of Cq and C. 

The grand partition function, evaluated in the saddle point approximation, can now be substituted in the standard 
thermodynamic relation in the presence of a chemical (electrostatic) potential 

S G =/3M-/3Q$+ \ogZ G , (14) 

so as to yield the grand canonical entropy 

Sq — Smc — 2 ^ cons t- (15) 

It is pertinent at this point to probe the validity of the saddle point approximation because it is used throughout 
this paper. One way to do this is to use this approximation to compute the mean-squared fluctuations of the energy 
and charge, 

AE 2 = (SE 2 ) = Z G X J dE dQ (E - M) 2 e SMC - pE+ ^ Q 

AQ 2 = (SQ 2 ) = Z G X J dE dQ (Q - Q ) 2 e SMC ~ (m+t ^ Q (16) 

and to ensure that these have certain desirable properties, e.g., that they are both positive and relatively small. 
Computing appropriate Gaussian integrals, it is easy to show that 3 

Smc,qq(M, Qq) 

det Q 
Smc,ee(M, Qq) 



AE 2 = 


- (fi- 


= 


AQ 2 = 


(JT 


X )22 = 


(6ESQ) = 


- (fi- 


"% = 



det n 

Smc,eq(M, Qq) 
det Q 



(17) 



Using the inequalities (11), (12) both mean-squared fluctuations are clearly seen to be positive, as they ought to 
be. The cross-correlation term, which is generally nonzero, explicitly demonstrates why it is erroneous to treat these 
fluctuations as independent, as already remarked. 

III. FLUCTUATION CONTRIBUTION TO THE ENTROPY OF LARGE BLACK HOLES 

The canonical black hole entropy has been argued [9] to be determined by the ADM Hamiltonian which describes 
quantum fluctuations on the horizon treated as an inner boundary of spacetime, essentially because the bulk Hamil- 
tonian obeys the quantum Hamiltonian constraint. Unfortunately, the spectrum of this boundary Hamiltonian is yet 



3 For a similar analysis of energy and particle number fluctuations in a grand canonical ensemble, see [15] 
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to be determined in LQG. Thus, as in [8], we assume that energy spectrum is a function of the discrete area spectrum 
(well-known in LQG [3] ) and a discrete charge spectrum. The charge spectrum is of course equally spaced in general; 
for large macroscopic black holes the area spectrum is equally spaced as well. 

In a basis in which both the area and charge operators are simultaneously diagonal, the grand canonical partition 
function can be expressed as 

Zg = Yl 9(m,n) exp-/3[ E(A m ,Q n ) - $ Q„] , (18) 

where, g(m,n) is the degeneracy corresponding to the area eigenvalue A m and charge eigenvalue Q n . Using a 
generalization of the Poisson resummation formula 



X] /( m > n ) = X] dx dy exp {- 

m.n k.l 



■i{kx + ly)} f(x,y) (19) 



and assuming that the partition sum is dominated by the large eigenvalues A m ,Q n , it can be expressed as a double 
integral 



Zg = J dxdy exp - (3 {E(A(x),Q(y)) - $ Q(y)} g 



(A(x),Q(y)). (20) 

Note that the transition from the discrete sum to the integral for Zg requires only that the dominant eigenvalues are 
large compared to the fundamental units of discreteness which for the area is the Planck area and for the charge is 
the electronic charge. These conditions are of course fulfilled for all astrophysical black holes. 
Changing variables in eq. (20) form x, y to E, Q 

Z G = I dE dQ J(E, Q) g(E, Q) exp {-0(E - $Q)} 

dE dQ p(E, Q) exp {-0(E - $Q)} , (21) 



where, the Jacobian J = \ E^ x | 1 | Q tV | 1 , and p = J(E,Q) g(E,Q) is the density of states. Employing the 
saddle point approximation and using eq. (14) one obtains 

S G {M,Q ) = Smc{M,Q ) - i log A + const. (22) 

where, using Smc = log/? = Smc + log J" we have defined 

A = dctftj 2 = detft (E, x ) 2 (g, y ) 2 |M,Q 

= dot ft (A B ) 2 (A iX ) 2 (Q, v ) 2 \m,q , (23) 

and £1 is defined in eq. (10). 

Since the microcanonical entropy is known to be only a function of the horizon area even for charged non-rotating 
black holes [4,5], one can express detO as 

detfl - {Smc.a? [(A.ee) (Aqq) - (Aeq) 2 ] \m,q 

+ Smc,aa Smc, a [{A,e) 2 A.qq + (A.q) 2 A^ee — 2 A^e A,q A,eq] \m,q (24) 

Now, to obtain the leading thermal fluctuation correction to the area law (over and above corrections occurring in 
Smc) it is adequate to retain only the area term in Smc in eq. (24), yielding 

dctft - ^ [AeeA, qq (A.eq) 2 } m ^ Qq ■ (25) 

With both the area and charge spectrum being equally spaced (i.e., linear in the quantum numbers x and y), we 
therefore obtain the leading thermal fluctuation correction 

A = ^ [A tEE A tQQ - (A.eq) 2 ] (A, E f |m,q (26) 
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where C is a constant related to the area and charge spectra. 

To proceed further, it is necessary to specify the functional dependence of the area of the horizon A on the energy 
E and the charge Q. As mentioned earlier, this is not yet known within the LQG framework. We adopt here the 
alternative of postulating an asymptotic power law relationship between these quantities for large A 

E = a A a + b Q 2 A~P , (27) 

where, a , b are real positive dimensional coefficients dependent on the Newton constant and the cosmological constant, 
and a , (3 are dimensionless real positive exponents. An equivalent way of rewriting eq. (27) is 



E = a A a 



1 + b'Q 2 A-?' , (28) 



where b' is similar to b above and likewise for j3'. Both forms of power-law relation subsume all four dimensional 
charged non-rotating black holes. Of these, only those spacetimes which are asymptotically anti-de Sitter, i,e, have 
constant negative curvature at infinity, can be considered to be in stable thermal equilibrium, so long as their (outer) 
horizon radius exceeds (— A) -1 / 2 where A is the negative cosmological constant. 4 

We may remark that there is no essential loss of generality in choosing the functional relation in eq.s (27, 28) to be 
quadratic in the charge; in fact our analysis goes through for any positive power of Q. 

Observe that the microcanonical ensemble definition of temperature is T = (Smc,e) ■ On physical grounds one 
expects that T > which implies a limit on the charge Q 2 < (aa/b/3) A a+ " using eq. (27). This inequality is 
saturated by the extremal charge Q 2 xt — (aa/b(3)A a+l3 for which T = 0. We also define, for future use, the ratio 
q = Qo/Qext which must satisfy — 1 < q < 1. 

The computation of the derivatives in eq. (26) is now straightforward, 

A,e\m,q = ^W(q,a,[3) A 1 -* , A, ee \ m .q = - ^ X(q,a,(3) A x ~ 2a 

Aqq\m,q = -\ Y(q,a,0) A 1 -"-" 
( b \ 1/2 

A,eq\m,q - f-H qZ(q,a,(3) A 1 - 3 ^- 1 ^ . (29) 

The dimensionless functions W, X, Y and Z are explicitly calculable and are regular and non- vanishing for < \q\ < 1; 
however they have poles at q = ±1 with W ~ (1 — q 2 )^ 1 and X , Y , Z ~ (1 — q 2 )~ 3 as \q\ — > 1. 
In terms of these functions, the leading correction is given by 

C b XY q 2 Z 2 (a+p) 

A - 16 a W 2 A ' (M) 

so that the grand canonical entropy 

S G - S MC + \ {pt + p) log A + 0(A°) for A » (-A)- 1 and < | 9 | < 1 . (31) 

Note that the leading logarithmic correction is independent of q. However, as we approach extremality \q\ — > 1, A 
diverges as (1 — q 2 )~ 4 , the O(A ) terms become large and the approximation of large horizon areas becomes invalid. 
This happens because as we approach extremality, T — > and all fluctuations basically freeze out. Furthermore, as 
we shall argue below, the result also cannot be trusted for q — > 0, although for somewhat different reasons. 

IV. RELATIVE FLUCTUATIONS FOR LARGE BLACK HOLES 

Using the area law, we find, for adS black holes with A » £ 2 , 



'in fact we only consider very large black holes with A S> (—A) for which (27) is assumed to be valid. 
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AE 2 = - 
AQ 2 = - 



1 (A, 



4 Vdct nj MiQo 

1 (Aee\ 

4 [det n) MQa 



(32) 



where, eq. (25) is to be used for detfi. Substituting the derivatives in eq. (29), we obtain the rms fluctuations in 
mass and charge 



AE 



AQ 




M,Q 



M — a A a ( 1 + —q 



Qa = q 



bp 



1/2 



On the other hand, we have for the mean (equilibrium) values 



and 



Hence the relative fluctuations are 

AE 

~w 

AQ 
Qo 



a 9 
1+ P q 



^("+/3)/2 



Y ) C-l/2 

XY - q 2 Z 2 J BH 



(33) 



(34) 



(35) 



(3X 



a(XY - q 2 Z 2 ) 



1/2 



? -l/2 
'BH 



(36) 



Thus, although the thermal fluctuation correction to the entropy has a degree of non-universality, the relative fluctu- 
ations nevertheless turn out fall off with A for large A in a universal manner, scaling like A~ x l 2 independent of the 
details of the energy-area relation. This universal fall-off shows that the small fluctuation approximation employed in 
our analysis is valid for large area black holes as long as q is not small. This requires the black hole charge to grow 
as a moderate proper fraction of the extremal charge which increases with area. 

However, for \q\ — > 0, while the relative rms energy fluctuation is bounded, the relative charge fluctuation blows up 
as q , thereby invalidating the small fluctuation (Gaussian) approximation for charge. This happens not only for 
the Qo ~ case considered in ref. [13] (where it also follows from the relation for the rms fluctuation in the charge 
quantum number, but appears to have been overlooked) but for any value of Qo which does not scale as Q ext (A) for 
large A. 



V. DISCUSSION 



The primary result of this paper is the leading thermal fluctuation correction ^ (a + (3) log A in the grand canonical 
entropy (31). It is obvious that when charge fluctuations are taken into account in addition to energy fluctuations of 
a black hole, the logarithmic correction due to such fluctuations depends on the exponents appearing in the relation 
linking horizon area, mass and charge. This is quite unlike the case studied in Ref. [8] where the charge fluctuations 
were ignored; the energy fluctuations in that case produced a correction \ log A independent of the assumed area- 
energy relation. In this sense, charge fluctuations endow the thermal fluctuation correction with a certain degree of 
'non-universality' absent in the earlier case. However, what is perhaps noteworthy is that the coefficient appearing 
in the leading thermal correction to the area law is proportional to the sum of the exponents appearing in the area- 
charge-energy relation (27), and not related to any one of the exponents alone. Another way of saying the same thing 
is that, if we use eq. (28) to express the area-charge-energy relation, the leading correction is in fact independent of a 
and simply depends on /?', which of course is the same as a + [3 . In (28), the contribution of the charge is expressed 
as an additional factor multiplying the energy-area relation used in the canonical ensemble case [8] . The fact that 
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the leading thermal correction is proportional only to the exponent appearing in this additional factor reaffirms our 
earlier result of universality in absence of charge fluctuations. To this extent then, one could say that the leading 
thermal correction with charge fluctuations is at least 'partly universal'. 

For the adS Reissner Nordstrom black hole, the area-charge-energy relation is given by 



where, I = (— A/3) -1 / 2 . For A » I 2 and moderate values of \q\ in the range < \q\ < 1 ), we have here a = 3/2 , f3' = 
2, in the notation used in (28, implying that the leading thermal correction is log A. This result agrees with that found 
in Ref. [13] for this black hole in the limit Qq — > (or equivalently q — > 0). This limit is however pathological since, 
as argued above, the mean-squared charge fluctuation relative to its mean (equilibrium) value diverges in this limit, 
thereby invalidating the Gaussian approximation which has been used to compute the thermal fluctuation correction. 
We have further demonstrated above that it is not at all necessary to go to this small charge limit in order to compute 
the leading thermal fluctuation correction. We may reiterate that our result eq. (31) actually encompasses a larger 
class of adS black holes which may include non-rotating adS dilatonic black holes as well [16]. 

It follows that the mechanism proposed in Ref. [14] to argue a possible complete cancellation of logarithmic cor- 
rections in the grand canonical computation of the entropy stands on loose ground. In contrast, such a cancellation 
has already been shown earlier [9] between the logarithmic contribution in the U(l) calculation of the microcanonical 
entropy [6] of non-rotating black holes and the universal canonical correction found by us earlier [8], both of which 
appear to be robust. Generalization of these considerations to black holes with rotation stands prominently on our 
agenda for the near future. 
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